Asymmetry in shape causing absolute negative mobility 
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We propose a simple classical concept of nanodevices working in an absolute negative mobility 
(ANM) regime: The minimal spatial asymmetry required for ANM to occur is embedded in the 
geometry of the transported particle, rather than in the channel design. This allows for a tremendous 
simplification of device engineering, thus paving the way towards practical implementations of ANM. 
Operating conditions and performance of our model device are investigated, both numerically and 
analytically. 
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I. INTRODUCTION 

Realizing a micro- or nano-device exhibiting absolute 
negative mobility (ANM) poses serious technological chal- 
lenges, as this task is believed to require finely tailored 
spatial asymmetries either in the (nonlinear) particle- 
particle interactions [l|, 0] or, more conveniently, in the 
geometry of the device itself 0, |3| . A device is said to op- 
erate in the ANM regime, when it works steadily against 
a biased force, i.e., a force with nonzero stationary mean. 
According to the Second Law of Thermodynamics (or 
more precisely, the so called principle of Le Chatelier), 
a static force alone cannot induce ANM in a device cou- 
pled to an equilibrium heat bath, unless an additional 
time dependent force is applied to bring the system out 
of equilibrium. ANM is known to occur as a genuine 
quantum mechanical phenomenon in photovoltaic mate- 
rials, as the result of photo-assisted tunneling in either 
the bulk of noncentrosymmetric crystals [|[ or artificial 
semiconductor structures 0, 0j- However, such manifes- 
tations of the ANM phenomenon do not survive in the 
limit of a classical description, so that detecting ANM in 
a purely classical system remains a challenging task. 

As spatial symmetry typically suppresses ANM, ad hoc 
contrived geometries have been proposed to circumvent 
this difficulty. The most promising solution devised to 
date, is represented by two (2D) or three dimensional 
(3D) channels with inner walls tailored so as to force the 
transported particles along meandering paths 0, H|, a 
design that can be implemented, e.g., in superconduct- 
ing vortex devices [8|. Other classical set-ups advocate 
elusive dynamic chaotic effects fij-fllj]. Although such 
finely tuned asymmetric geometries and/or nonlinear dy- 
namic behaviors may seem hardly accessible to table-top 
experiments, first convincing demonstrations of classical 
ANM have actually been obtained following this strategy 

mm. 

We propose here a much simpler, affordable working 
concept for a classical ANM device, by embedding the 
spatial asymmetry into the shape of the transported par- 
ticles, rather than in the channel geometry. In view of 
this new formulation, the ANM mechanism is expected 



to occur in natural systems, too, where cylindrically sym- 
metric channels in low spatial dimensions and elongated 
particles are frequently encountered [H], HU • 

This paper is organized as follows. We introduce in 
Sec. |TI]the Langevin equations for a floating ellipsoidal 
Brownian particle ac-driven along a 2D compartmental- 
ized channel. By numerical simulation we show in Sec. 
IIIII that ANM actually occurs as an effect of the parti- 
cle elongation. In Sec. IIVI we analyze the dependence of 
ANM on both the drive parameters and the particle ge- 
ometry, with the purpose of determining the optimal op- 
erating conditions of our model device. Finally, in Sec. |V] 
we discuss the applicability of the proposed ANM mecha- 
nism to nanoparticle transport in realistic biological and 
artificial devices. 
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FIG. 1: (Color online) (a) Asymmetric particle tumbling in 
a periodically segmented 2D-channel. The pores, 2A wide, 
are centered on the channel axis, (b) Elliptic particle with 
semiaxes a and 6, at rest against a compartment wall. Note 
that a < A < b. (c) Elliptic particle in escape position, its 
major axis forming a maximum angle (f>* with the channel 
axis. 
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II. THE MODEL 

Let us consider an elongated Brownian particle, shaped 
as an elliptic disk, moving in a straight 2D channel (Fig. 
[T]) . The overdamped dynamics of the particle is modeled 
by three Langevin equations, namely 

— = -F(t) e x + JW r ((t) ; (la) 

^ = y/Dit*(t), (lb) 

where e x ,e y are the unit vectors along the x, y axes, 
r = (x, y) denotes the particle center of mass, and 4> is the 
orientation of its major axis with respect to the channel 
axis, e x . Here, £(f) = (^ x {t)^ y (t)) and ^(t) are zero- 
mean, white Gaussian noises with autocorrelation func- 
tions (&(*)&(*')) = 25ij5(t - t') and i,j = x,y,<p. The 
channel is periodically segmented by means of orthogonal 
compartment walls, each bearing an opening, or pore, of 
half- width A, placed at its center As sketched in 

Fig. QJa), the channel is mirror symmetric with respect 
to both its longitudinal axis and each compartment wall. 
This is an important difference with Ref. [H, HJ , where 
the channel confining potential, V(x, y), was taken to 
be asymmetric under both mirror reflections - although 
symmetric under double reflection, V(—x, y) = V(x, — y). 

In order to detect ANM, the particle must be driven 
in a pulsating manner parallel to the channel axis. This 
means that F(t) consists of at least two terms [1,3: a dc 
drive, Fq, and an unbiased, symmetric ac drive, F ac (t), 
with amplitude m&x{\F ac (t)\} = Fi and temporal period 
Tq. Accordingly, the waveform of F ac (t) is subjected to 

the conditions {-Fic" (t))n = 0, with n = 0, 1, 2 . . . and 
(. . . )q denoting the time average taken over one drive 
cycle (l7j . 

Equations ([TJ have been numerically integrated for an 
elliptic disk of semiaxes a and 6, under the assumption 
that the channel walls were perfectly reflecting and the 
particle- wall collisions were elastic [lg]. In the following 
we report the outcome of extensive simulations for a fixed 
channel compartment geometry, xl = Vl, ^lyh *C 1, 
but different particle elongations, b/a, ac drive wave- 
forms, and ratios of the rotational to translational dif- 
fusion coefficients, D^/D,.. We conclude that ANM oc- 
curs in such a highly symmetric channel geometry only 
because of the elongated aspect ratio of the drifting par- 
ticle. 

As illustrated in panels (b) and (c) of Fig. [I] an elliptic 
disk with a < A < b crosses a narrow pore only when 
its major semiaxis forms a small angle with the channel 
axis, \(f>\ < (j>*. For a <C b (rod- like particle) as in most of 
our simulations, sin</>* ~ A/6. To overcome the escape 
angle (f>* from a rest position with <p — 7r/2, the disk must 
rotate against the total applied force, F(t) = F + F ac (t). 
For Fq < Fi this is more easily achieved for pore crossings 
occurring in the direction of F ac (t), but opposite to the 
static force Fq. As a result, under appropriate conditions, 
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FIG. 2: (Color online) ANM for a driven-pulsated elongated 
Brownian particle: current (v) vs. static bias Fq in the pres- 
ence of square-wave drives with amplitude F\ = 2 and dif- 
ferent periods Th (see legend). Other simulation parameters: 
Diffusion strengths D — D r — D$ =0.1; shape parameters 
a = 0.05, b — 0.3; compartment parameters xl = ])l = 1, and 
A = 0.1. Each data point for (v) = linii^oo (x(t) — x(0))/t 
was computed from a single trajectory with t = 10 6 and time- 
step 10 -5 ; the statistical error was estimated to be 5%, i.e., 
of the order of the symbol size. Note, for a comparison, that 
the compartment traversal time is tq — 7 and the diffusive 
relaxation times are = = 5. The dotted curve repre- 
sents v(F) at zero ac-drive, F\ = 0. Inset: the corresponding 
mobility curve (solid curve), /J.{F), is compared with the an- 
alytical estimate of Eq. ((4]) (dashed curve) for F m = 0.86. 



detailed below, the net particle current, (v), may indeed 
flow in the direction opposite to Fq. 



III. THE ANM MECHANISM 

In order to explain the appearance of ANM, we start 
looking at the mobility of an elongated particle driven by 
a constant force F (Fig. ^ inset). Let v(F) denote its 
steady velocity and fi(F) = v(F)/F the relevant mobil- 
ity, with n(—F) — fJ<(F). From now on, and until stated 
otherwise, we set for simplicity D r = = D, so that 
fi is a function of F/D. At equilibrium with F = 0, 
fio = /i(0) is a Z?-independent constant, which strongly 
depends on both the compartment and the particle ge- 
ometry as discussed below. For zero drive, the particle 
rotates away from the walls; pore crossing is thus con- 
trolled mostly by translational diffusion. 

As the magnitude of the applied force is increased, the 
mobility of elliptic and circular disks develop a quite dif- 
ferent F dependence. The mobility of a circular disk 
with a = b < A, n(F), is a concave function of F/D, 
which decays from fig = (x(0) to fj,^ = fi(F — > oo) = 
(A — a)/(yL — a), with a power law slower than 
[l8|]. In the case of an elliptic disk with a < A < 6, 
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reaching the escape angle <f>* can be regarded as a noise 
activated process with energy barrier proportional to F. 
Pore crossing will then be controlled mostly by the rota- 
tional fluctuations, with approximate escape time 

T (F)=T exp(F/F m ), (2) 

where F m = 2D/ (b cos <f>* — a) is the total activation force, 
the factor 2 accounts for the two directions of rotation, 
and To is the compartment traversal time, xl/F, divided 
by the probability, p = (A — cl)/(jjl — a), that the disk 
slides through the pore without an additional rotation. 
The reciprocal of r plays the role of an effective attack 
frequency. This estimate for the particle crossing time 
surely holds good for b » A and F 3> F m , where 

^o«^ ) ,r^ ) «T (F), (3) 

with = x 2 L /2D and = y 2 L /2D denoting, respec- 
tively, the longitudinal and transverse relaxation times. 

The curve /i(F) for an elongated particle is thus con- 
cave for F < D/xl [13 and decays exponentially for 
F » F m , like 

H(F) - = 2pexp (-F/F m ), (4) 

see inset in Fig. ® Correspondingly, v(F) increases like 
HoF at small F and decays to zero like xl/t(F) at large 
F, going through a maximum for F ~ F m , as conhrmed 
by our simulations, see Fig. [2] 

A convincing evidence for ANM has been obtained by 
simultaneously applying to the elliptic disk a tunable 
dc force, Fo, and a low- frequency, square- wave ac force, 
Fac(t)t with amplitude F\ > Fo > 0. The characteristics 
curves (v) vs. F plotted in Fig. [5] exhibit a negative 
ANM branch only for sufficiently long ac drive periods 
and F < Fx; for F > Fx, however, the current is always 
oriented in the Fq direction, no matter what Tq. 

This behavior can be explained in the adiabatic regime, 
where a half drive period, To/2, is larger than all the 
drift and diffusion times inside a channel compartment, 
namely, r , , and . Note that in the opposite 

regime, ANM is suppressed. The net current can then be 
approximated by 

(v(Fo)) = \[v{Fx + F ) - v(Fx - F )]. (5) 

As the curve v(F) peaks at F ~ F m , we expect (v(F )) 
to develop a negative minimum for Fq = F\ — F m and a 
positive maximum for Fq = Fx + F m , as shown in Fig. [21 
As this holds true only for Fx > F m , the two peaks have 
upper bounds 

\(v(Fx - F m ))\ < {v(Fx + F m )) < v(F m ) . (6) 

For a more quantitative analysis of this phenomenon, 
we rewrite v(Fx±F ) = (Fx ± F ) fi(Fx ± F ) , so that the 
ANM condition, (v) < 0, reads 

F l > + F o) „ p -2F /F m m 

F 1+ F > »(Fx-F )- e ■ {<) 




FIG. 3: (Color online) (a) Drive waveform dependence: (v) 
vs. Fo for three waveforms of F ac (t), square (as in Fig. [2]), 
sinusoidal and ramped, all with F\ = 2, D — 0.02, and Tn = 
10 3 . (b) Particle inertia dependence: (v) vs. Fo for three 
values of m = I, see text. Other simulation parameters are 
Fx = 2, D = 0.1, and T n = 500. 

The approximate equality on the r.h.s. applies for < 
F < Fx - F m , where 

M(F 1 ±F )~ M (F 1 )e^o/^. (8) 

The above inequality is satisfied for < Fo < F* , with 
the turning point, F* , shifting towards zero in the limit 
Fx — > F m +, and towards Fx in the opposite limit, Fi 3> 
F m . The approximate equality in Eq. ([7]) leads to slightly 
overestimating F*, with no prejudice of our conclusion: 
In the adiabatic regime, ANM occurs in an appropriate 
F interval (0,F*) only provided that Fx > F m . 

Finally, we notice that for Fo > Fx > F m the net 
current reads 

(v(Fo)) - \[v{F + Fx) + v(F - Fx)} (9) 

and for extremely large F , it decays to zero like 
(v(F )) ~ ^v(F — Fx). Correspondingly, our simulation 
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FIG. 4: (Color online) Particle elongation dependence: {v) 
vs. b for different values of Tn (main panel) and of D (in- 
set). Other simulation parameters: square- wave ac-drive with 
amplitude strength F\ = 2, and bias Fo = 1, a = 0.05, 
xl = y.L = 1, and A = 0.1. 



data in the neighborhood of the turning point F* ~ Fi 
are reasonably well reproduced by the linear fitting law, 
M^o)}~(/V2)(Fi -*!>). 



IV. SELECTIVITY AND OPTIMIZATION 
CRITERIA 



In view of future experimental implementations of the 
proposed ANM mechanism, we now analyze in detail its 
sensitivity with respect to both the drive and the particle 
parameters. 

We start noticing that the most prominent ANM ef- 
fect is produced, in fact, by the square waveform F ac (t) 
adopted in Fig. [2J For the sake of a comparison, in 
Fig. EJa) we plotted (v(Fq)) also for other, inversion- 
symmetric waveforms F ac (t) with the same amplitude 
and period, in particular, sinusoidal and up-down ramped 
waveforms. For a ramped ac drive, no ANM can occur, 
because in the adiabatic regime 



I rF 1 +F 

M*b)> = 2ljr / v(F)dF>0 

1 J Fi — Fq 



(10) 



For a sinusoidal ac drive, the ANM effect can be shown 
analytically to diminish in magnitude and shrink to a 
narrower interval (0,F*) than obtained for the corre- 
sponding square waveform. The latter is thus the optimal 
ac drive waveform to operate an ANM device. 

To quantify the robustness of this effect against the 
damping conditions, inertia was added to the model by 
replacing the l.h.s. in the Langevin equations Q as fol- 
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~di 
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In Fig. |3]Jb) we compare ANM characteristics curves for 
growing values of the (rescaled) particle mass, m, and 
moment of inertia, /: ANM is gradually suppressed by 
increasing inertia. This is no serious limitation, as in 
most experiments rectifiers operate, indeed, under over- 
damped, or zero mass, conditions [8|. 

We analyze next how selective the ANM effect is versus 
the geometric and diffusive properties of the transported 
particles. In Fig. 2] we displayed the dependence of the 
net current on the particle elongation. One notices im- 
mediately that, when plotted versus b at constant values 
of the drive parameters, (v) starts out positive and then 
turns negative for b larger than a certain threshold, b*, 
which appears to increase with either raising D (figure 
inset) or lowering Tq (main panel). 

Our adiabatic argument provides a simple explanation 
for these findings, as well. We recall that the mobility 
curve /J,(F) decays exponentially on a scale F m oc b^ 1 . 
As a consequence, for b — > oo, Eq. (O boils down to 
(v) ~ — \v(F\ — Fq), which means that (v) tends to zero 
from negative values, in agreement with our data. For 
b < A the particle flows through the pores, no matter 
what the orientation, 4>, of its major axis. Therefore, 
(v) becomes insensitive to the particle elongation (main 
panel), while retaining its known D dependence (inset). 

The actual value of b* is determined by the general 
ANM condition ([7]). On making use of the approxima- 
tion on the r.h.s. of that equation, one easily proves the 
existence of the threshold b* for any geometry and drive 
parameter set. We caution that this way one may un- 
derestimate b* and, therefore, the predicted dependence 
b* oc D holds qualitatively, only (see inset of Fig. |4}. Of 
course, when D is raised so that b* grows larger than the 
spatial dimensions of a channel compartment, then ANM 
is suppressed altogether. 

The dependence of the current on the fluctuation in- 
tensities is illustrated in Fig. [5] We consider first the case 
D r = D t jj = D (figure inset). The dependence of (v) on 
D can be analyzed following the approach introduced to 
interpret the results in Fig. 2) On recalling that F m oc D, 
in the limit D — > oo, the ac drive amplitude ends up being 
smaller than F m , F m > Fi, thus suppressing ANM. In 
the opposite limit, D — > 0, F m vanishes and ANM is pre- 
dicted to occur for any dc drive such that < Fo < Fi. 
Indeed, from Eq. ([5]) we obtain [i(Fi ± Fq) — > fi or 
(v) — > (iqFo > 0, for D — s> oo (marked in figure by 
horizontal arrows), and (v) ~ — \v{F\ — Fo) — > 0—, for 
D — > 0. On using D as a control parameter, ANM is thus 
restricted to low noise, < D < D* , with the threshold 
D* also obtainable from the ANM condition (J7J. 

We consider next the more general case when D r and 
Dtj, can be independently varied, while keeping a and b 
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FIG. 5: (Color online) Rotational translational diffusion de- 
pendence: (v) vs. D^jD r for different values of D r (see leg- 
end). Inset: (v) vs. D = D r = with D$/D r = 1. The 
arrows mark the asymptotes predicted in the text. Other 
simulation parameters are as in Fig. [4] with b — 0.3. 

fixed. In Fig. [3J (v) has been plotted versus D^jD r 
for different values of D r . The two opposite limits of 
the net current, (v)o, for D<p/D r — >• 0, and (w)oo, for 
D^/ D r — >• oo, are both positive with (v)o < (w)oo. In be- 
tween, the magnitude of the ANM effect is seemingly not 
much sensitive to D^,/D r over several orders of magni- 
tude, which allows us to generalize the conclusions drawn 
above for D^/D r = 1 to the case of realistic extended 
particles. Note that (v)q is non-null because an elliptic 
disk with a < A < b can diffuse across a pore even in 
the limit — > 0, where the adiabatic argument fails, 
thanks to the sole translational fluctuations. In view of 
the crossing condition \<p\ < </>*, an elongated particle 
can be handled as a circular one with radius smaller that 
A, see Fig. 0] but crossing probability 2<p* /tt. This ar- 
gument can be extended to the case of (w)oo, with the 
important difference that for — > oo the particle has 
crossing probability one. Both (v)o and {v)oo are thus 
positive, with (v}q relatively smaller than (w)oo- 

V. CONCLUDING REMARKS 

The ANM model presented in this work, although styl- 
ized, lends itself to interesting nano-technological appli- 
cations [1, HH . A typical reference case is represented by 
the transport of hydrated DNA fragments across narrow 
compartmentalized channels [20| . Rod- like DNA frag- 



ments 30-40nm in length and l-2nm in (hydrated) di- 
ameter are easily accessible [21| ; their elongation ratio is 
about 3 times larger than b/a in Fig. [2j but still within 
the ANM range of Fig. 01 Artificial nanopores can be 
TEM drilled in lOnm thin Si02 membranes with repro- 
ducible diameters of 5nm, or less [201 > which is consistent 
with the elongation selectivity condition a < A < b as- 
sumed throughout this work. Moreover, the measured 
Dd/Dr ratio for the hydrated DNA fragments of Ref. 
21] falls in the range 100-200, in dimensionless units, 
where ANM can also occur, as shown in Fig. for an 
appropriate choice of the drive parameters. To this re- 
gard, we remind that experiments on DNA translocation 
across artificial nanopores require applied electrical fields 
of the order of 10-100kV/cm 22]; if applied to the DNA 
rods of Ref. [2l|, electrical fields of that intensity, or 
less, would satisfy the ANM condition of Eq. ([7]) with 
F\ > F rn at room temperature. 

The ANM characteristics curves plotted in Figs. [2] and 
[5J being quite selective with respect to the particle shape, 
suggest the possibility of developing artificial devices that 
efficiently operate as geometric sieves for nanoparticles. 

Our model was stylized to capture the key mecha- 
nism responsible for the occurrence of ANM in symmetric 
channels. The mechanism summarized by Eqs. ([2]) and 
([7]), however, clearly does not depend on the dimension- 
ality of the channel (experiments can then be carries out 
in 3D geometries), but can be impacted by other compet- 
ing effects: (i) Pore selectivity. For a given translocating 
molecule, the actual crossing time varies with the wall 
structure inside the pore and in the vicinity of its open- 
ing [23j]; (ii) Electrophoretic effects. The inhomogeneous 
electrical field generated by the electrolyte flow across the 
pores acts on the orientation of drifting spheroidal parti- 
cles [24| • System specific effects (i) and (ii) can readily be 
incorporated in our model by adding appropriate poten- 
tial terms, U r (x, y) and U,p(4>), to the Langevin equations 
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